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We present a numerical study of the quantum action previously introduced as a parametrisation 
of Q.M. transition amplitudes. We address the qüestions: Is the quantum action possibly an exact 
parametrisation in the whole range of transition times (0 < T < 00)? Is the presence of potential 
terms beyond those occuring in the classical potential required? What is the error of the parametri- 
sation estimated from the numerical fit? How about convergence and stability of the fitting method 
(dependence on grid points, resolution, initial conditions, infernal precision etc.)? Further we com- 
pare two methods of numerical determination of the quantum action: (i) global fit of the Q.M. 
transition amplitudes and (ii) flow equation. As model we consider the inverse square potential, for 
which the Q.M. transition amplitudes are analytically known. We find that the relative error of the 
parametrisation starts from zero at T = increases to about 10 -3 at T = 1/E gr and then decreases 
to zero when T — > 00. Second, we observe stability of the quantum action under variation of the 
control parameters. Finally, the flow equation method works well in the regime of large T giving 
stable results under variation of initial data and consistent with the global fit method. 
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I. WHAT IS THE QUANTUM ACTION? 



In Refs.|l|]-|5| the quantum action has been introduced. It concerns the foundations of quantum physics. 
It basically states that transition amplitudes in quantum mechanics can be expressed in terms of an action, 
which has the form of the classical action, but with parameters (mass, potential) which are different from 
those of the classical action. This presents a new link between quantum physics and classical physics. 
The quantum action can be interpreted as a renormalized action in Q.M. |5|. The quantum action is 
similar to the Standard effective action, however, it is free of the deseases of the latter (infinite series of 
higher derivative terms, non-localities etc). In the limit of large imaginary transition time, the existence 
of the quantum action has been proven [Q. In this limit, the quantum action has a number of remar kable 
properties: (i) The WKB approximation f3CÍ|| for the ground state wave function becomes exact, after 
replacing the classical action by the quantum action. (ii) There is a differential equation relating the 
classical mass and potential to the quantum mass and potential. (iii) There is an analytic expression for 
the ground state wave function in terms of the quantum action. (iv) The ground state energy coincides 
with the minimum of the quantum potential. (v) The ground state wave function has a maximum exactly 
at the same position, where the quantum potential has a minimum. (vi) The quantum action allows also 
to reproduce energies and wave functions of excited states. Example: The spectrum of the hydrogen atom, 
considering the lowest energy states for given angular momentum. 
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FIG. 1: Inverse square potential. Classical potential with parameters m = 1, uj = 1 and g = 1(T 2 (full line), 
quantum potential (dashed line) and ground state wave function (enhanced by factor 5, dotted line). 



2 



Inverse square potential 



6.5 
6 

5.5 
5 

0.508 
0.506 
0.504 
0.502 
0.5 
1 

0.995 
0.99 
0.985 



Singular term 



Harmònic term 



Mass 



FIG. 2: Inverse square potential. Parameters of the quantum action rh (bottom), V2 (center) and v~2 (top) as 
function of transition time T. Classical parameters m = 1, 1)2 = 0.5 and V-2 = 5. 



II. USE OF THE QUANTUM ACTION 

Like the Standard effective action || , the optimized expansion of effective action M , the Gaussian effective 
action [||, the Feynman-Kleinert effective action M also the quantum action has been constructed for the 
purpose to study quantum phenomena, which have its origin in classical physics. Prototype examples are 
chaos and instantons. ít is well known that classical chaos has no direct analogue in quantum physics. One 
reason is that in Q.M. one cannot define a point in phase space (due to Heisenberg's uncertainty relation), 
hence one cannot define Lyapunov exponents from diverging trajectòries, thus quantum chaos cannot be 
represented quantitatively via Poincaré sections. ín Ref.Q the effective action has been used to describe 
classically chaotic quantum systems. A comparison of Poincaré sections from the classical action and the 
quantum action has been presented in Ref. q{ . A problem similar to that of a proper definition of quantum 
chaos exists also for the definition of instantons in quantum physics. Classical instantons are solutions 
going from one potential minimum to another. They start from the location of the potential minimum 
and with kinetic energy zero. Due to Heisenberg's uncertainty relation, there is no rigorous analogon in 
quantum physics. The physics of instantons and its relation to tunneling has been discussed in Refs.pïj |ï^| . 
For a review on the use of the Standard effective action see Ref. [131^ A comparison of instantons from the 
classical action versus the quantum action has been given in Ref. 0] . Because tunneling in Q.M. is very 
closely related to instantons, the quantum action is expected to shed new light into the phenomenon of 
Q.M. tunneling. Below we give a survey of tòpics in different areas of physics where the use of the different 
kinds of effective action as well as the quantum action should be useful. 



A. High energy physics 



(i) Renormalisation. The quantum action may turn out to provide a new definition of renormalisa- 
tion. In conventional terminology, renormalisation mcans to extract physical observable parameters like 
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FIG. 3: Inverse square potential. Parameter m of quantum action vs. transition time T. Quantum action is 
parametrized by supplementary terms x 4 and x~ 4 . Classical parameters m = 1, vv = 0.5 (lo — 1), w_2 = 1 (g = 1). 
Boundary points of transition: xí - 2 points in interval [4, 5]; Xf - 10 points in interval [0.5, 3]. Temporal resolution 
Aí = 2 10" 3 . 



mass, coupling etc. obtained from an interacting quantum field theory at the continuum limit. Those 
physical parameters difFcr from the so-called bare parameters and tiris difference represents the effects of 
the interaction. Very similar to this, the quantum action has parameters (mass, potential parameters) 
diffcrcnt from the classical action and this difference represents the effects of quantum mechanics: Q.M. 
fluctuations oceuring in the path integral are represented by a single path, however, for a particle with 
different properties of mass and interaction. In order to explore the use of the quantum action as means of 
renormalisation, one could construct the quantum action for a many-body system (like a chain of coupled 
anharmonic oscillators) and explore its quantum mechanical continuum limit. In quantum field theory, 
one carries out renormalisation by computing n-point vèrtex functions, which represent vacuum-to-vacuum 
transition amplitudes for transition time T — > oo. Exactly in the limit T — > oo, the existence of the 
quantum action has been proven rigorously in Q.M., although so far only in the case of a single particle 
system. 

(ii) Cosmology and inflationary scenario of early universe. Inflation involves potentials with several 
minima and instanton solutions. The instanton starts out as a quantum instanton and eventually turns 
into a classical instanton. This has effects on the subsequent formation of galàxies 0, |Ï5|, [ïq| . Using 
the effective potential or the quantum potential (potential of the quantum action) in general creates a 
potential different from the classical one. In particular, it may have minima being absent in the classical 
potential. Consequently, this may create instantons being quite different from the classical instanton (in 
Ref. H instantons from the quantum action were found to be "softer" than the corresponding classical 
instanton). Such quantum effects of the instanton may influence the outeome of the galaxy formation at 
the end of inflation. 

(iii) Hot and dense nuclear matter. Instantons are believed to play an important role in hot nuclear 
matter in the quark-gluon plasma phase |Ï7|]. Instantons are important also for the mechanism of chiral 
symmetry breaking and for t'Hoofts solution of the U(l) problem. 

(iv) Neutrino oscillations. The process of oscillations of neutrino flavors |ï^, |Ï9L ^Ü) may have to do with 
a process of tunneling in a potential with degenerate minima. 
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FIG. 4: Same as Fig.Q for parameter 
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FIG. 5: Same as Fig.Q for parameter 



V-2. 



B. Condensed matter physics 



(i) Quantum dots, semi-conductor and quantum chaos. Advancing the speed of microprocessors may 
have technological obstacles but also physical límits. When reducing the size of a chip one soon may enter 
the regime where quantum laws rule. Quantum chaos may become a very important issue, because it 
can hamper the flow of elèctric currents. The quantitative determination of quantum chaos effects will 
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FIG. 6: Same as Fig.H for parameter v- 




FIG. 7: Same as Fig.Q for parameter 114. 



possibly be of great importance for the development of future microprocessors. By use of some kind of 
effective action or the quantum action, one can study quantum corrals formed by àtoms and quantum dots 
in semi-conductors. In particular, this allows to study the temperature dependence of electron dynamics 
in atòmic corrals, as well as for electrons moving in simple conductor-semiconductor-isolator geometries 
and to search for the possible presence of quantum chaos. 

(ii) Josephson junctions and superconducting quantum interference. Superconducting quantum inter- 
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FIG. 8: Same as Fig.H for parameter V-4. 
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FIG. 9: Same as Fig.j3|. Global relative error of quantum action Ey. 



ference devices (SQUID) have been used to demonstrate experimentally the phenomenon of quantum 
superposition in macroscopic states [pï| . This involves Josephson junctions. The SQUID potential has a 
double-well structure. The effective potential and the quantum potential should be useful tools to analyze 
quantum superposition in terms of such potential involving quantum effects. 

(iii) Quantum computers based on superconductors. The symmetry of the order parameter in some 
triplet superconductors corrersponds to doubly degenerate chiral states. Gulian et al. 123, 03] predict that 
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FIG. 10: Inverse square potential. Mass rh vs. transition time T. Dependence of parameters of quantum action on 
location of final boundary points: Xf - 100 points in various intervals. Initial boundary point Xi = 0.3. 
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FIG. 11: Same as Fig.]ÏÒ]| for parameter üq. 

this degeneracy can be lifted via macroscopic quantum tunneling. Instanton-like quantum behavior may 
become important. Triplet superconducters may be used as bàsic elements of quantum computers. Again 
the effective potential and the quantum potential should help to study quantum instantons and tunneling 
in such materials. 
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FIG. 12: Same as Fig.jï^ for parameter V2- 
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FIG. 13: Same as Fig.JÏ^] for parameter C_2. 
C. Atòmic physics 

(i) Analogy of classical chaos in quantum physics. Attempts have been made to use the effective action in 
order to characterize chaotic behavior in quantum systems [ÏC| ] . In a a similar way, the quantum action has 
been used also H . This allows to construct of a phase space portrait and Poincaré sections for a quantum 
system in analogy to classical physics. From this one can obtain Lyapunov exponents and KAM surfaces 
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for the quantum system. A potentially quite interesting system to explore is the Paul trap or similar traps. 

(ii) Ultracold àtoms in a billard formed by làsers. Trace formulas (Gutzwiller |Q and generalisations ) 
have been used successfully [^5|, |27j to establish a relation between level densities and periòdic semi- 
classical orbits. It would be interesting to compare predictions of trace formulas in the semi-classical regime 
with the predictions obtained from the effective action or the quantum action. 

(iii) Dynamical tunneling. Steck et al. (28) and Hensinger et al.[^9| have demonstrated experimentally 
the phenomenon of dynamical tunneling (where the classical transition is forbidden due to some conserved 
quantity different from energy). It has been realized by arrays of cold àtoms. It has been observed that 
the presence of quantum chaos enhances the dynamical tunneling transition. It would be instructive to 
reexamine dynamical tunneling using the phase space portrait constructed from a time-dependent effective 
action or quantum action. 

D. Chemistry 

Binding of macromolecules. In the process of chemical binding of macromolecules, often a double well 
potential plays a role. The effective action as well as the quantum action should be useful to find pathways 
in the formation of such macromolecules. 



III. QUANTUM ACTION VS. EFFECTIVE ACTION 

What are the similarities of the Standard effective action and the quantum action? Where are the 
differerences? 

(i) The effective action is computed via loop expansion. This leads to non-localities. The effective 
action is given by an infinite series of terms. Those non-localities are manifested by higher order terms of 
time-derivatives [ (^)™, (^f )™\ etc] The series of such higher terms is infinite. This causes problems 
because the solution of the Euler-Lagrange equation of motion requires infinitely many boundary conditions. 
Moreover, such series does not converge, it is only asymptiotically vàlid. To make practical use of the 
effective action one proceeds by retaining a few low order terms. However, for the purpose to study 
quantum chaos |ïc(| , this may be problemàtic, because in chaotic dynamics "small terms" may cause large 
cffccts. None of such non-localities oceur in the quantum action. 

(ii) While the Standard effective action does not reproduce a potential of double well shape, the Gaussian 
effective action produces an effective potential with double well shape |?], É|. This is also the case for the 
Feynman-Kleinert type of effective action [|| . 

(iii) Q.M. tunneling has no classical counterpart. It is closely related to instantons. One should note 
that the particular tunneling process and also the instantons depend on the shape of the potential, its 
barrier height as well as the position of minima. Because different types of effective action differ in their 
corresponding effective potential, consequently they differ in their tunneling amplitudes. 

(iv) The quantum action shares with the effective action the property that the ground state energy 
is given by the minimum of the quantum, respectively effective potential. But the quantum action is 
distinguished by the property of coincidence of location(s) of maximum(a) of the ground state wave function 
and minimum(a) of the quantum potential j5j. 

IV. VALIDITY OF THE QUANTUM ACTION 

The quantum action has been originally proposed as a conjecture. A priori, it is not evident that such 
quantum action exists. By now the existence of the quantum action has been established in the following 
cases: 

(i) Harmònic oscillator. In this case the quantum action is identical to the classical action. 
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FIG. 14: Inverse square potential. Mass rh vs. transition time T. Dependence of parameters of quantum action on 
location of initial boundary point Xi. Final boundary points Xf - 100 points in interval [2, 3]. 
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FIG. 15: Same as Fig.Q for 
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(ii) In the limit when the transition time T — > 0, the quantum action exists. Dirac has noticed that 
the path integral in this limit is dominated by the contribution from the classical trajectory. Hence the 
quantum action coincides with the classical action. 

(iii) In imaginary time (necessary to describe thermodynamics at finite temperature) and going to the 
limit T — > oo, the quantum action exists. In general it is quitc diíferent from the classical action. 
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FIG. 16: Same as Fïg.@ for 
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FIG. 17: Same as Fïg.Q for í- 



This leads to the question: What about finite transition time? Does the quantum action exist for 
arbitray finite transition time? Does it parametrize Q.M. transition amplitudes for all Xi n to Xfi equally 
well? In general, the quantum action has to be determined non-perturbatively by numerical computations. 
What can be said about about stability, convergence and errors of such procedure? Those qüestions are 
the subject of the paper. We have chosen to study this at hand of the inverse square potential in one 
dimension. 
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FIG. 18: Inverse square potential. Mass m vs. transition time T. Dependence of parameters of the quantum action 
on boundary points: Xi - 10 points in interval [1.5, 2.5]; Xf - 10 points in interval [1.1, 2.1]. 



V. INVERSE SQUARE POTENTIAL 



We consider the following classical potential in 1-D 

V(x) — —muj 2 x 2 + g x~ 2 



v 2 x 



V-2 x 



(D 



The corresponding classical action is given by S — J dt ^mx 2 — V(x), while the classical action in imaginary 
time (í — > —it) is given by Se = J dt ^mx 2 + V(x) (note: following conventions used in physics, the 
overall minus sign has been dropped. It reappears in the weight factor exp[— Se])- The potential is parity 
symmetric. Because it has an infinite barrier (for g > 0) at the origin, the system at x < is separated 
from the system at x > 0. We consider only the motion in the domain x > 0. The potential is shown in 
Fig.[Q. We have chosen to consider the inverse square potential, because of the distinct feature that the 
corresponding quantum mechanical transition amplitudes are known analytically |33|, |34| . They are given 
by 



G{b,T;a,0) = 
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where L· is the modified Bessel function and 
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The numerical studies discussed below have been done in imaginary time. The Euclidean transition am- 
plitude reads 



G E (b,T;a,0) 



ab 



hsmh(ujT) 



exp 



ftsinh(u;T) 



(4) 
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FIG. 19: Same as Fig.[|18| for «o- One observes vo — > 0.5, consistent with analytic result (minimum of quantum 
potential goes to E gr ). 
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FIG. 20: Same as Fig.Jïit for mV2- One observes mvi — * 0.5, consistent with analytic result. 



The transition amplitude contains all information on the spectrum and wave functions. For example, by 
going to the limit T — > oo (Feynman-Kac limit), the Euclidean transition amplitude is projected onto the 
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FIG. 21: Same as Fig. |lq] for mü-2. One observes that m»_2 — > 2, consistant with analytic result. 
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FIG. 22: Same as Fig.@. Global relative error of Gi 



ground state. Onc finds 



G E (b, T; a, 0) -^ T ^oo Z exp[-w(l + f)T\ 



,1/24-y r mUJ 

b L/2+1 exp 1 
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l/2+i r mw 2i 



(5) 
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where 

2muj fmujy 

Zo = jmTT){—) ■ (6) 

One reads off the ground state energy and wave function 

E gr — hid(l+-f) 

= 4 /2 * 1/2+7 ^M-^A ■ (7) 

This is, of course, identical with the direct solution from the Schròdinger equation. The wave function is 
shown in Fig.[Q. 
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FIG. 23: Inverse square potential. Mass rh vs. transition time T. Dependence on the density of meshpoints N per 
unit interval of T. 



A. Dynamical time and length scales 

In this work we do a numerical study of a model to test the quantum action. Our choice of the model 
has been infhicnccd by the analytical solvability and less by the question if it plays a role in nature. 
Consequently, absolute vàlues in physical units of the model parameters (mass, potential parameters) are 
not of primary importance. We have expressed those parameters in dimcnsionless units. However, what 
is important are time and length scales, which are dynamically generated by the model, say a time scale 
T sc and a length scale A sc . Those scales serve as reference vàlues to give a sense to statements like "for 
large transition times T", or a "small spatial resolution Ax" , which means T/T sc >> 1 and Ax/A sc << 1, 
respectively. How to choose those dynamical scales? We have used as time scale 
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FIG. 24: Same as Fig-fea] for v Q . 
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FIG. 25: Same as Fig.feí for v 2 . 



As a length scale one may introduce the analogue of the Bohr radius of the ground state wave function. 
Another possibility is to define a length scale A sc by 



dx \ip gr (x)\ 2 = 0.95 



(9) 
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FIG. 26: Same as Fig.kÈ for w_ 



i.e. the length which covers 95% of the probability of the ground state wave function. For the classical 
action considered in Eq. po| of sect.[Vl], those scale parameters are 



T sc = 0.4 
A sc w 2.35 



(10) 



B. Analytical results for the quantum action in the asymptotic regime 

In Ref.g we have shown for the Euclidean asymptotic regime T — > oo, that the following analytic 
relations exist between ground state energy E gr , ground state wave function ip gr and the quantum action. 



(11) 



i.e. the minimum of the quantum potential gives the ground state energy. Next, the ground state wave 
function can be expressed in terms of the parameters of the quantum action 



- _1 r-ÍL·- dX ' ^fl™(.V(x')-V mm )/h 

N 



(12) 



Finally, there is a relation between classical and quantum action (transformation law) which reads (for 



2m(V(x) - E gr ) = 2m(V(x) - V rnm ) -'- dx 



H £2m(V(x) - V min ) 



2rh(V(x) - V min ) 



(13) 



Here x m i n denotes the position of the minimum of the quantum potential. Let us make an ansatz for the 
quantum action, characterized by a mass fa and a quantum potential of the form 



V(x) = V2X 2 + V-2X 2 + v 



(14) 
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and see if it satisfies transformation law, Eq.(p"3|). Using Eqs.PJÏ4|), inscrting it into Eq.([Ï3|), and comparing 
the coefncients of the terms x 2 , x~ 2 and x°, one obtains 



2m V2 = m lü 



2 , ,2 



Those equations arc equivalent to 



2m ü_2 — fr\/2rh v-2 — 2m g 
4m\Jv2 v-2 + hy / 2m V2 — 2m E gr 



rh V2 — 2 171,2 w 

m w_ 2 = ^h 2 [~ + 7 ] 2 

B ffr = n [i + 7] ■ 



(15) 



(16) 



For the case when the classical potential is given by the parameters m=l,/ï=l,cj = l and g = 5, this 
yields 



m u_2 
rh í>2 — 



-+ 0.5 
6.8507 



The numerical solution, shown in Fig.[2] gives 



m ü_2 -> 0.500 
rh V2 —* 6.86 . 



(17) 



(18) 



We observe that analytical and numerical results agree well. 

Finally, let us compare the analytic behavior of the classical function mV(x) (mass x potential) with the 
corresponding quantum function rhV(x). In particular, let us look at the term x~~ 2 which is singular at the 
origin. The corresponding term in the classical potential is mg/x 2 (Eq.(0). Let us see what happens when 
g — > 0. The classical term mg/x 2 is singular as a function of x for any g > 0. For g = the singular term 
vanishes and the classical potential is regular as a function of x. This means that the classical potential 
has a singular behavior as a function of g when g — > 0. Now let us look at the corresponding term in the 
quantum potential Eq.([Ï6|) telis us 



m V-2 



= h 2 
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2 n 2 
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Amg 



0(g 2 ) 



(19) 



This function has a regular behavior when g — > 0. Thus, in the limit g — > 0, the quantum potential 
is smoother than the classical potential. A similar behavior has been observed previously in a study 
comparing a classical double well potential with the corresponding quantum double well potential and its 
corresponding instanton solutions g] . The quantum double well potential was found to have a lower barrier 
and the quantum instantons were found to be softer than their classical counterparts. 



VI. NUMERICAL RESULTS FROM GLOBAL FIT 



Let us first consider the classical action with the following parameters 

m = 1 , V2 = 0.5 , v_2 = 1 • 

This gives the ground state energy 

Eg r = 2.5 . 



(20) 
(21) 
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For the quantum action the following ansatz has been madc, 



S = I dt -mx 2 + V(x) 

V(x) — Vq + V2X 2 + V-%X~ 2 + V4X 4 + V-4X 



(22) 



The method to determine the parameters of the quantum action from given transitions amplitudes 
G(b, T;a,0) proceeds by selecting a number of initial positions x^ and also final positions x$ . Mak- 
ing a guess for the parameters of the quantum action, one computes the trajectory corresponding to the 
quantum action with the boundary condition x(t = 0) = x^ and x(t = T) = x^). Then one computes 
the value of the quantum action for this trajectory, say T, pa . This is done for all o~ and p. The goal is to 
make the error of the fit 



globally small. The parameters of the quantum action are those which minimize the error. For more details 
on the method see Ref. [Q. In the calculations reported here we have used the following initial and final 
data: There are 2 initial points Xi n located in the interval [4, 5]. There are 10 final points Xfi equidistantly 
distributed in the interval [0.5,3]. For the computation of the trajectòries, the following resolution of time 
has been used: N t = 500 meshpoints for a time interval AT = 1. For larger time intervals, the number 
of meshpoints has been increased proportionally. An overview of the behavior of the parameters of the 
quantum action is displayed in Figs.|}]- D- Fig.p shows the mass of the quantum action m as a function 
of transition time, which displays a smooth behavior. The behavior of the estimated errors seems to be 
smooth between T = up to T = 4. At T = 4, the error bars behave erratically. The behavior of the 
parameters vq, V2 etc. also looks smooth. In this calculation we included terms x and x in the quantum 
potential. The data show that the coeffients are quite small and almost compatible with the value zero. 
Fig.|5) presents the relative error of the quantum action vs. transition time T. The error starts out at 
T = from a value close to zero, then increases and reaches at T = 1.5 a maximum of 5 10 -4 and finally 
decreases until T = 5. At T = 5 there is a cusp. By adding terms x~ 4 and x 4 to the quantum potential 
gives a qualitatively similar behavior. It reduces the maximum of the error by a factor of about 2. Note 
that the parameters are shown on a linear scale while the error of the quantum action is shown on a 
logarithmic scale in order to resolve smaller numbers. In order to see if higher order terms play a role in 
the quantum potential, we have carried out a calculation including the terms x 6 and x~ 6 instead of x 4 and 
x~ 4 . The results (not shown here) are quite similar to the previous ones. 

Dependence on boundary points. Next we look at the dependence of the quantum action parameters on 
the boundary points used in the fit. We have studied this in two ways: (a) We kept the initial points Xi 
fixed and varied the final points Xf. (b) We kept the final points Xf fixed and varied the initial points 2^. 
In the first case, we have considered one initial point Xi n = 0.3 (kept fixed). As final points Xfi we took 
100 points uniformly distributed in an interval and varied those intervals [2,3], [5,6], [9,10] and [2,10]. 
The results are shown in Figs.0- One observes for all parameters of the quantum action that in the 
regime T > 2 (T > 5T SC ) there is no dependence on the location of the final boundary points. However, 
for T < 2 (T < 5T SC ) there is a noticeable dependence. 

In the second case, as final boundary points xt we took 100 points uniformly distributed in the interval 
[2, 3] (interval kept fixed). As initial points xi we took one point and varied it between 0.1 to 0.5. The results 
are shown in Figs.jÏ4| - |l7| ]. Those results are qualitatively the same as in the first case. Fig.jÏ4| shows 
that the quantum mass depends on the boundary points for T < 2 (T < 5T SC ) and becomes independent 
for T > 2 (T > 5T SC ). The behavior for the parameters of the quantum potential is similar, however, the 
independence sets in a bit earlier at about T sa 1.6 = 4T SC . 

One may ask: What is the reason for such dependence on boundary points? In our opinion the most 
likely explanation is that the quantum action is not an exact parametrisation of the transition amplitude 
in this regime. Hence the parametrisation depends on the parameters of the fit. 




(23) 
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The two previous cases represent a strong imbalance between the number of initial versus final boundary 
points. As a third case we have considered a more balanced selection of bondary points. We took 10 
initial points, uniformly distributed in the interval [1.5,2.5] and 10 final points, uniformly distributed in 
the interval [1.1, 2.1]. The parameters of the quantum action are shown in Figs.|ÏJ-|2Ï). The global relative 
error in fitting the transition matrix elements is shown in Fig.j2^]. 

Dependence on temporal resolution. The determination of the parameters of the quantum action by a 
global best fit, Eq.(|2j), requires to solve the equation of motion from S for all pairs of boundary points 
and to compute the value S of the quantum action along such trajectory. Such numerical calculations 
depend on the temporal resolution Aí. We have studied the convergence of the parameters of the quantum 
action as a function of the density of meshpoints. We parametrize the density of meshpoints by taking 
N t meshpoints per unit time interval AT = 1. We have varied N t from 200 up to 8000. The results are 
shown in Figs.j2^ - Fig.j23) shows for T < 6 » 15T SC that a low density Nt = 200 is sufficient to reach 
convergence. However, for larger T, the results diverge. If one desires convergence for, say T = 8 ~ 20T SC , 
one needs to double the density of meshpoints. Such behavior persists if we want to maintain convergence 
for even larger T. E.g., for T = 14 w 35T SC , we need a 20-fold higher density of meshpoints (Nt — 8000). 
Such exponential increase in the density of meshpoints could signal chaos caused by imperfect numerical 
solutions, fmitc infernal prccision, rounding errors etc. 

Dependence on internat precision. We have looked at the parameters of the quantum action in the 
asymptotic regime (T — > oo) and attempted to analyze to role of infernal computing precision. We 
performed a computation using FORTRAN giving approximately 20 significant digits and a computation 
using MAPLE giving approximately 30 significant digits. We found (not shown here) at T w 15 « 38Ts C 
a bifurcation between the computation with 20 and 30 digits. Moreover, the computation with 20 digits 
becomes flat for T > 15, while the computation with 30 digits becomes flat for T > 19. Thus, a tiny 
change in the infernal precision by e — 10 -30 results in a change of m at T = 19 by a margin of 2.5%. This 
signals a great sensibility to infernal precision. Again, this could be interpreted as the presence of chaos, 
due to mathematical algoritms. (One should note that the TD Hamiltonian system under consideration is 
integrable and hence physical chaos is absent). 



VII. RENORMALISATION GROUP EQUATION FOR TEMPERATURE DEPENDENCE OF 

ACTION PARAMETERS 

A global fit method has been used in sectfVÏ| to determine the parameters of the quantum action. The 
method is quite costly from the computational point of view. It requires first to compute the quantum 
transition matrix elements and then to carry out a multi-parameter fit. Here we want to suggest an alterna- 
tive non-perturbative method. It is inspired by the idea of the renormalisation group equation in quantum 
field theory. In QFT the dependence of action parameters upon variation of a scale parameter (cut-off 
A, lattice spacing a s and at) is governed by a renormalisation group equation (e.g. Callan-Symanzik). In 
the Q.M. system considered here, we are close to the continuum limit (in the example considered below 
Ax/ob = 0.045 and At/as = 5 x 10 -6 , where as is the ground state Bohr radius of the ground state wave 
function). 

Here we draw a parallel between the scale-dependence of an observable action parameter in field theory 
(e.g. coupling g in QCD) and the dependence of the parameters of the quantum action upon (Euclidean) 
transition time T. Finite temperature physics is obtained from transition amplitudes in Euclidean transition 
time T, and the inverse temperature (3 is related to T via (3 — T/h. As scale and temperature dependence of 
the action are similar, we apply here the term renormalisation group to describe temperature dependence 
of the quantum action ||. As a result we will end up with a differential equation (flow equation) for 
the parameters of the quantum action. In contrast to the global fit method, this flow equation does not 
require as input the transition amplitudes. It does, however, require initial data, i.e. the parameters of the 
quantum action at some initial temperature. 
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FIG. 27: Quantum action parameters for quartic potential. Comparison of results from flow equation (lines) vs. fit 
to transition amplitudes (symbols). 



We consider the transition amplitude as a function of x and t, keeping initial data Xi n , ti n fixed. It 
satisfies the Schròdinger equation 



H G(x, t, Xi n , ím) 



2m dx 2 



V{x) 



G(x,~t] Ximtin) i 



with the initial condition 



lim G(x,t;x in ,ti n ) = S(x - x in ) 



(24) 



(25) 



In the limit t — > í in , the transition amplitude is given by the classical action, consistent with the initial 
condition, Eq.(||), 



t lim ^G(x,t;x in ,t in ) - ^exp[-^Sc/|^ jti j| -> 



(26) 



Going over to inverse temperature /3, the parametrisation of (Euclidean) transition amplitudes in terms of 
the (Euclidean) quantum action reads 



G = Zp expj-Èp] 



(27) 
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Combining Eq.(J27|) with Eq.(|24D (expressed in terms of (3) implies 

1 dZ fj dT, a h 2 r ,dtfi^ dPt 



-0 



Z dP df3 



fi 2 r/ d£, 

— ( — - 

2m LV dx 



Y 



dx 2 



V = 



(28) 



Up is given by the quantum action along its classical trajectory from Xi n , f3i n = to x, (3 



£tf = S, 



0[X c l 



ix,/3 



x,/3 



,dx 



tp [m(p),vo(p),vi(p) t . 



2h 2 K d/3 
,x,0_ 



(29) 



Here Vk(P), k = 0, 1, 2, . . . denotes a set of parameters of the quantum potential (as an example Vk((3) 
may be the coeficients of the quantum potential V{x,0) expanded in terms of polynomials) . In general 
the number of terms is infinite. The weight of higher terms decreases rapidly. Eqs . (p9| , p8| ) yicld the 
renormalisation group equation for the action parameters, m(/3), vq(/3), v\{(3), ... as function of j3, 



1 dZp dY.fi d 



' 



Z Q d/3 drh 8(3 



E 



OY, 



dv k d(3 d(3 



h 2 

— [( 

2m dx 



dYp^ 2 d 2 T, l3 ^_ v = Q 



?- 



dx 2 



(30) 



This equation is vàlid for all x, Xi n . The parameters rh(/3), Vk((3) are independent of x, Xi n . This constitutes 
a system of equations to determine dZp/d(3, drh/ d (3 and dvk/d/3. The solution of the differential equations 
requires initial vàlues. We noted the property that the quantum action goes over to the classical action 
in the limit when the transition time goes to zero. This suggests to take for [3 = as initial vàlues 
—> 0) ~ Z(0 ~ * 0) (note: singularity at origin), rh(fi = 0) = m, Vk((3 = 0) = Vk, k — 0,1, ... . 
Starting from initial vàlues the renormalisation group equation determines the flow of the renormalized 
parameters when (3 increases. 
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FIG. 28: Inverse square potential. Quantum action parameters obtained from flow equation. Dependence on initial 
data of flow equation. Mass m vs. transition time T. 
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Inverse square potential 




FIG. 29: Same as Fig.M for v Q . 
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FIG. 30: Same as Fig.pJ for v 2 . 



Flow equation for quartic potential. As a first test of the method we have applicd it to the quartic 
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FIG. 31: Same as Fig.M for iL 



potential in 1-D, given by the classical action 



S 



dt —mx 2 
2 



V(x) 



V{x) = V2X + V4X 

m = 1, 1)2 = 1, v 4 = 0.01 



(31) 



We have computed numerically the solution of the flow equation of the quantum action parameters, Eq. (pü|) . 
We have taken a single initial boundary point Xi n and considered transitions to a number of final boundary 
points Xf. We have compared this with a global fit of the quantum action to QM transition amplitudes 
(see Ref.jj]). A comparison is shown in Fig.|27). One observes good agreement in the range (5 — 0.5 to 
= 4. 

Flow equation for inverse square potential. We considered the action with the inverse square potential, 
with the parameters of the classical action given by 



1, ,v 2 — 0.5, u_2 = 1 



(32) 



which are the same parameters as used in Eq.(J20|) of sect.(Vl). We looked at the behavior of the parameters 
of the quantum action obtained by solving the flow equation. First, we want to study the stability of the 
parameters of the quantum action under variation of initial conditions. In particular we chose the following 
initial conditions: At Ti n u — 0.35, we used as initial data parameters obtained by the global fit method, 



m 



0.99994533, v = 1.1676274, £_ 2 = 1.2280919 



(33) 



The inital value of V2 has been varied between 0.47 to 0.52. In particular, we have studied how the quantum 
action parameters depend on the choice of initial value of «2- As initial boundary points we used a single 
point, Xi n — 10. As final boundary points Xf we used 30 points uniformly distributed in the interval [0.2, 7]. 
The differential equations have been solved with a resolution A/3 — 3.75 10 -3 . The results are shown in 
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Figs.|p8| - ^ï[ . One observes for rh and í>_2, shown in Figs.[^8| |3~ï| |, a very weak influence of the variation 
of the initial value of V2- However, in vq, Fig. [p9| , there is a substantial variation. It is interesting to note 
that V2, shown in Fig. |30j, converges rapidly for (3 > 2. In other words, except for vq, all parameters of the 
quantum action, although starting from different initial data, collapse to a single curve, when (3 > 2. The 
flow equation method shows stability of results under variation of initial data. 

Second, we looked at the dependence of the solution on the location of final boundary points. In this 
case, we took at Ti„n — 0.35 the following initial data 

m = 0.99994533, v = 1.1676300, v 2 = 0.4998810, w_ 2 = 1.2280900 . (34) 

We took a single initial boundary point xi n = 10, As final boundary points Xf we took 30 points uniformly 
distributed in some interval [c, çfL where c = 0.2 and d has been varied from d = 5 to d = 17. The 
results are displayed in Figs.|Ï2| - B3j. Fig.|3^] shows a weak dependence on the location of final boundary 
points. There is some visible dependence for T < 2. For T > 2, this dependence is invisible from the 
figure. However, because we have an analytic prediction for rnü_2 in the limit t oo, we can compare 
the numerical results with the analytic prediction. The difference is plotted in Fig. |33]. One observes that 
some dependendence on the location of the final boundary points continues to exists beyond T = 2 (note 
logarithmic scale). However, in the limit of large T, this difference tends to zero (within numerical errors). 
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FIG. 32: Inverse square potential. Quantum action parameters obtained from flow equation. Dependence on 
location of interval [c, d] of final boundary points. rhv-2 vs. T. 



VIII. DISCUSSION 

We have studied a 1-D Hamiltonian system, which is integrable, given by an inverse square potential. It 
has the attractive feature that its quantum transition amplitudes are analytically known. We have chosen 
this system to test the validity of the quantum action. The knowledge of transition amplitudes eliminates 
one source of error in the construction of the quantum action. We carried out a numerical study and 
determined the parameters of the quantum action by two different methods: via global fit and via flow 
equation. 
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We consider the following observations as most important: As a function of transition time T (or inverse 
temperature 0), the behavior of the parameters of the quantum action can be divided in three regimes: (i) 
Regime of small T, say < T < T sc . For small T high lying energies will contribute in thermodynamic 
observables. In atòmic physics, the high-lying states of the hydrogen atom are known as Rydberg states 
and the wave functions are known to have a semi-classical behavior. Thus this regime is the semi-classical 
regime. This is consistent also with the well known property that quantum transition amplitudes in the 
limit T — ► are well described by G — Z exp[jrS c i], i.c. the path integral of the transition amplitude is well 
approximated by the classical action evaluated along a single path - the trajectory of the classical action. 
It is interesting to note that trace formulas (Gutzwiller and extensions) work quite well in the semiclassical 
regime. Then the path integral can be well approximated by periòdic orbits. In this work we observe that 
the quantum action also works well in the semiclassical regime. More precisely, it works the better the 
smaller the value of T is. 

(ii) There is another regime, that of large T, say 5T SC < T < oo. This regime is opposite of the 
semiclassical regime, i.e. the deep quantum regime. In thermodynamics, large (3 means small temperature. 
In this limit the Feynman-Kac formula holds and the physics is described by the ground state properties. 
It is not known if the trace formulas work in this regime (the authors are not aware of any evidence) . It is 
known that the quantum action becomes an exact parametrisation in this limit . Here we have analyzed 
the behavior of the quantum action numerically. We found that both, the global fit method, as well as the 
flow equation method give results in agreement with each other and also with the analytic result. 

(iii) Finally, there is an intermediate regime, say T sc < T < 5T SC . This is the regime where the 
parametrisation of the transition amplitude gives the largest relative global errors (see Figs.f^, ||), also 
manifested by the strongest dependence on boundary points. Although this error is small, it is by several 
orders of magnitude larger than the error in the regime of large T . How can this be understood? The 
following scenarios are possible: First, the ansatz for the quantum potential may be incomplete. There may 
be a need to include local terms in the quantum potential beyond those occuring in the classical potential. 
We have carried out some steps in this direction. We included in the quantum potential the following 
terms, all absent from the classical potential: x 4 and x~ 4 . Some improvement has been found (see Fig.|^]), 
however, it is less than one order of magnitude. Similar results have been obtained when incorporating the 
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terms x e and x~ 6 . Of course it is possible that there are terms missing of a forní quite different from the 
polynomial form considcrcd here. However, we consider that as unlikely. The results seem to indicate that 
this error is not mainly due to a lack of terms in the quantum potential. 

Second, the parametrisation of transition amplitudes by the quantum action may be incomplete in the 
sense that it is not sufhcient to evaluate the quantum action along its corresponding classical trajectory, 
which minimizes the quantum action, but other trajectòries need to be taken into account. Let us be 
more specific. For stiff potentials, like the quartic potentials, it is well known that there exists an infinity 
of classical trajectòries, for a given pair of boundary points. This has been discussed in some detail 
by Schulman |3Ï], |3^]. Those trajectòries are all solutions of the classical equations of motion, and all 
correspond to the same pair of boundary points, but they differ in the value of its action. On the other 
hand, we observed that the quantum action works well in the limit of large T, i.e. the Feynman-Kac regime, 
where the physics is given by the ground state properties. When lowering T then gradually the first excited 
state, the second excited state, etc. will contribute in the partition function and the (Euclidean) transition 
amplitudes. At the same time we notice an increase of the error in the parametrisation by the quantum 
action when lowering T . We propose the following hypothes: In a regime where higher lying states need 
to be taken into account in the transition amplitudes, then the parametrisation by the quantum action 
nccds to take into account trajectòries of the quantum action with higher vàlues of the action. There 
is an interesting geometrical analogy between the wave functions of excited states and those "higher" 
trajectòries: The ground state wave function has no nodes, but the wave functions of excited states have 
nodes (its number increasing with energy). Correspondingly "higher" trajectòries have wiggles (its number 
increases with the value of the action). 

Third, the concept of the quantum action as a means to parametrize transition amplitudes may be plain 
wrong. But then one has to ask: Why does the quantum action work in the regime of small and large T? 
In our opinion the most promising route to further explore those qüestions is to carry out high precision 
numerical simulations in different models. 
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